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^ Abstract 

> : 

Recent developments in the study of primordial black holes (PBHs) are reviewed, with partic- 
[ ular emphasis on their formation and evaporation. It is still not clear whether PBHs formed but, if 
they did, they could provide a unique probe of the early Universe, gravitational collapse, high energy 
| physics and quantum gravity. Indeed their study may place interesting constraints on the physics 
>• ' relevant to these areas even if they never existed. 

m : 

t^. , Preface 

' It is a great honour to speak at this meeting because Katsuhiko Sato has been a friend and 

O colleague for 25 years. I first met him when he visited Cambridge in 1980. At the time he was a 
i- £h ^ ! guest professor at NORDITA and I was a research fellow at the Institute of Astronomy. It was a 
i particularly important period for Katsu because he had just written his famous paper on inflation 
^ (arguably the first paper on the topic, although its publication was delayed). I remember our first 
^ 1 meeting very vividly: he came to my office and asked to take my photograph because he said I was 
1 famous. Nobody had ever told me I was famous before, so I was quite chuffed! Of course, Katsu is 
now much more famous than me, so I intend to repay the compliment by taking a photograph of him 
?h ' at this meeting! 

Although we have not written any papers together, Katsu was the first Japanese physicist with 
whom I interacted and he has had an important influence on my career. This is because meeting him 
was one of the factors which prompted me to spend a year at Kyoto University in 1982 and I have 
made many subsequent visits to Japan. On several of these visits I have enjoyed the hospitality of 
Katsu and his wife, so I would like to express my gratitude for this. Indeed I would say that Katsu's 
scientific contribution lies not only in his own work but also in the role he has played in encouraging 
collaboration with Western scientists (many of whom are attending this conference). In my case, this 
collaboration is particularly strong, since two years ago I married a Japanese lady. It is strange to 
reflect that this might not have happened without the chain of events triggered by Katsu coming to 
my office in 1980! 

I've chosen to talk about primordial black holes, partly because this is one of my own lifelong 
interests, but also because it is an area in which Japanese physicists have made a vital contribution. 
Katsu has himself written many papers on the subject |132l 11531 11041 1105} 11271 19*T1 11541 1155j and I 
will highlight some of these in this paper. PBHs also bear on inflation, the topic of the meeting: on 
the one hand, the fluctuations which arise at inflation are the most likely source of PBH; on the other 
hand, upper limits on the number of PBHs place constraints on inflationary models. 
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1. Introduction 

Black holes with a wide range of masses could have formed in the early Universe as a result 
of the great compression associated with the Big Bang [7H1 1171j . A comparison of the cosmological 
density at a time t after the Big Bang with the density associated with a black hole of mass M shows 
that PBHs would have of order the particle horizon mass at their formation epoch: 



PBHs could thus span an enormous mass range: those formed at the Planck time (10 _43 s) would have 
the Planck mass (10 _5 g), whereas those formed at 1 s would be as large as 1O 5 M , comparable to the 
mass of the holes thought to reside in galactic nuclei. By contrast, black holes forming at the present 
epoch could never be smaller than about 1M Q . 

The realization that PBHs might be small prompted Hawking to study their quantum prop- 
erties. This led to his famous discovery [ZH] that black holes radiate thermally with a temperature 



Only black holes smaller than about 10 15 g would have evaporated by the present epoch, so eqn (1) 
implies that this effect could be important only for black holes which formed before 10~ 23 s. 

Hawking's result was a tremendous conceptual advance, since it linked three previously dis- 
parate areas of physics - quantum theory, general relativity and thermodynamics. Even if PBHs never 
existed, it has been useful to think about them! However, at first sight it was bad news for PBH 
enthusiasts. For since PBHs with a mass of 10 15 g would be producing photons with energy of order 
100 MeV at the present epoch, the observational limit on the 7-ray background intensity at 100 MeV 
immediately implied that their density could not exceed 10 -8 times the critical density |142j . Not 
only did this render PBHs unlikely dark matter candidates, it also implied that there was little chance 
of detecting black hole explosions at the present epoch |146j . Nevertheless, it was soon realized that 
the 7-ray background results did not preclude PBHs playing other important cosmological roles (2S] 
and some of these will be discussed here. 

2. How PBHs Form 

The high density of the early Universe is a necessary but not sufficient condition for PBH 
formation. One also needs density fluctuations, so that overdense regions can eventually stop ex- 
panding and recollapse. Indeed one reason for studying PBH formation and evaporation is that it 
imposes important constraints on primordial inhomogeneities. PBHs may also form at various phase 
transitions expected to occur in the early Universe. In some of these one require pre-existing density 
fluctuations, but in others the PBHs form spontaneously, even if the Universe starts off perfectly 
smooth. In the latter case, limits on the number of PBHs constrain any parameters associated with 
the phase transition. The formation of PBHs from inhomogeneities is discussed in Section 4, but we 
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Soft equation of state. Some phase transitions can lead to the equation of state becoming 
soft for a while. For example, the pressure may be reduced if the Universe's mass is ever channelled 
into particles which are massive enough to be non-relativistic. In such cases, the effect of pressure 
in stopping collapse is unimportant and the probability of PBH formation just depends upon the 
fraction of regions which are sufficiently spherical to undergo collapse jHE] • For a given spectrum of 
primordial fluctuations, this means that there may just be a narrow mass range - associated with the 
period of the soft equation of state - in which the PBHs form. 

Collapse of cosmic loops. In the cosmic string scenario, one expects some strings to self- 
intersect and form cosmic loops. A typical loop will be larger than its Schwarzschild radius by the 
factor (G/i) -1 , where \x is the string mass per unit length. If strings play a role in generating large-scale 
structure, G\i must be of order 10~ 6 . However, as discussed by many authors [1451 l%T| IBT^ 1231 . 
there is always a small probability that a cosmic loop will get into a configuration in which every 
dimension lies within its Schwarzschild radius. This probability depends upon both [i and the string 
correlation scale. Note that the holes form with equal probability at every epoch, so they should have 
an extended mass spectrum. 

Bubble collisions. Bubbles of broken symmetry might arise at any spontaneously broken sym- 
metry epoch and various people have suggested that PBHs could form as a result of bubble collisions 
jUJ |H3 I105| 11131 H33j . However, this happens only if the bubble formation rate per Hubble volume 
is finely tuned: if it is much larger than the Hubble rate, the entire Universe undergoes the phase 
transition immediately and there is not time to form black holes; if it is much less than the Hubble 
rate, the bubbles are very rare and never collide. The holes should have a mass of order the horizon 
mass at the phase transition, so PBHs forming at the GUT epoch would have a mass of 10 3 g, those 
forming at the electroweak unification epoch would have a mass of 10 28 g, and those forming at the 
QCD (quark-hadron) phase transition would have mass of around 1M . There could also be worm- 
hole production at a 1st order phase transition |1U4[ I127j . Recently the production of PBHs from 
bubble collisons at the end of 1st order inflation has been studied extensively by Khlopov and his 
colleagues |TTU1ITUU] . 

Collapse of domain walls The collapse of sufficiently large closed domain walls produced at a 
2nd order phase transition in the vacuum state of a scalar field, such as might be associated with 
inflation, could lead to PBH formation jlHJ I148j . These PBHs would have a small mass for a thermal 
phase transition with the usual equilibrium conditions. However, they could be much larger if one 
invoked a non-equlibrium scenario |149j . Indeed Khlopov et al. argue that they could span a wide 
range of masses, with a fractal structure of smaller PBHs clustered around larger ones |101j . 

3. Why PBHs are Useful 

The study of PBHs provides a unique probe of four areas of physics: the early Universe; 
quantum gravity; gravitational collapse and high energy physics. One can probe the last two topics 
only if PBHs exist today but one can gain insight into the first two topics even if PBHs never formed. 

PBHs as a probe of the early Universe (M < 10 l5 g). Many processes in the early Universe could 
be modified by PBH evaporations |14()j . For example, recent studies have considered how evaporations 
could change the details of baryosynthesis jlH] and nucleosynthesis |132l HU6j , provide a source of 
gravitinos [HH] and neutrinos swallow monopoles [9"Tlll59j and remove domain walls by puncturing 
them [lfiOj . PBHs evaporating at later times could have important astrophysical effects, such as 
helping to reionize the universe [H21- There has also been interest in whether PBHs could preserve 
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gravitational "memory" if the value of G was different in the early Universes [U ^1 EHJ E21 EH] • 
Other constraints are associated with thermodynamics |115j and the holographic principle [4"o] . 

PBHs as a probe of gravitational collapse (M > lO 15 ^. This relates to recent developments in 
the study of "critical phenomena" and the issue of whether PBHs are viable dark matter candidates. In 
the latter case, their gravitational effects could lead to distinctive dynamical, lensing and gravitational- 
wave signatures. They could also influence the development of large-scale structure and the formation 
of supermassive black holes in galactic nuclei. These issues are discussed in Section 5. 

PBHS as a probe of high energy physics (M ~ 10 15 g ). PBH evaporating today could contribute 
to cosmic rays, whose energy distribution would then give significant information about the high 
energy physics involved in the final explosive phase of black hole evaporation. In particular, PBHs 
could contribute to the cosmological and Galactic 7-ray backgrounds, the antiprotons and positrons in 
cosmic rays, gamma-ray bursts, and the annihilation-line radiation coming from centre of the Galaxy. 
These issues are discussed in Section 6. 

PBHS as a probe of quantum gravity (M ~ lO -5 ^. Many new factors could come into play 
when a black hole's mass gets down to the Planck regime, including the effects of extra dimensions 
and quantum-gravitational spacetime fluctuations. For example, it has been suggested that black 
hole evaporation could cease at this point, in which case Planck relics could contribute to the dark 
matter. More radically, it is possible that quantum gravity effects could appear at the TeV scale 
and this leads to the intriguing possibility that small black holes could be generated in accelerators 
experiments or cosmic ray events. Although such black holes are not technically "primordial", this 
would have radical implications for PBHs themselves. These issues are discussed in Section 7. 

4. Formation of PBHs from Inhomogeneities 

One of the most important reasons for studying PBHs is that it enables one to place limits on 
the spectrum of density fluctuations in the early Universe. This is because, if the PBHs form directly 
from density perturbations, the fraction of regions undergoing collapse at any epoch is determined 
by the root-mean-square amplitude e of the fluctuations entering the horizon at that epoch and the 
equation of state p = ■yp (0 < 7 < 1). One usually expects a radiation equation of state (7 = 1/3) in 
the early Universe but it may have deviated from this in some periods. 

4-1. Simplistic Analysis 

Early calculations assumed that the overdense region which evolves to a PBH is spherically 
symmetric and part of a closed Friedmann model. In order to collapse against the pressure, such a 
region must be larger than the Jeans length at maximum expansion and this is just y/j times the 
horizon size. On the other hand, it cannot be larger than the horizon size, else it would form a 
separate closed universe and not be part of our Universe |2Zj • 

This has two important implications. First, PBHs forming at time t after the Big Bang should 
have of order the horizon mass given by eqn (1). Second, for a region destined to collapse to a 
PBH, one requires the fractional overdensity at the horizon epoch 5 to exceed 7. Providing the 
density fluctuations have a Gaussian distribution and are spherically symmetric, one can infer that 
the fraction of regions of mass M which collapse is j21] 

(4) 



P(M) ~ e(M) exp 
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where e(M) is the value of e when the horizon mass is M. The PBHs can have an extended mass 
spectrum only if the fluctuations are scale-invariant (i.e. with e independent of M). In this case, the 
PBH mass distribution is given by [23] 



where M* « 10 g is the current lower cut-off in the mass spectrum due to evaporations, VLpbh is the 
total density of the PBHs in units of the critical density (which depends on 3) and the exponent a is 
determined by the equation of state: 



a = 5/2 if one has a radiation equation of state. This means that the density of PBHs larger than 
M falls off as M _1//2 , so most of the PBH density is contained in the smallest ones. 

Many scenarios for the cosmological density fluctuations predict that e is at least approxi- 
mately scale-invariant but the sensitive dependence of 8 on e means that even tiny deviations from 
scale-invariance can be important. If e(M) decreases with increasing M, then the spectrum falls off 
exponentially and most of the PBH density is contained in the smallest ones. If e(M) increases with 
increasing M, the spectrum rises exponentially and - if PBHs were to form at all - they could only 
do so at large scales. However, the microwave background anisotropies would then be larger than 
observed, so this is unlikely. It is also possible to enhance PBH formation at particular mass-scales 
by introducing bumps in e(M). 

The current density parameter f2pBH associated with PBHs which form at a redshift z or time 
t is related to 8 bv [24] 



where Qr ~ 10 -4 is the density parameter of the microwave background and we have used eqn (1). 
The (1 + z) factor arises because the radiation density scales as (1 + z) A , whereas the PBH density 
scales as (1 + z) 3 . Any limit on I^pbh therefore places a constraint on 3(M) and the constraints 
are summarized in Fig. El which is taken from Carr et al. [321 • The constraint for non-evaporating 
mass ranges above 10 15 g comes from requiring fipBH < 1 but stronger constraints are associated with 
PBHs smaller than this since they would have evaporated by now. The strongest one is the 7-ray limit 
associated with the 10 15 g PBHs evaporating at the present epoch |142j . Other ones are associated 
with the generation of entropy and modifications to the cosmological production of light elements 
|140j . The constraints below 10 6 g are based on the (uncertain) assumption that evaporating PBHs 
leave stable Planck mass relics, an issue which is discussed in detail in Section 7.1. 

The constraints on 3(M) can be converted into constraints on e(M) using eqn (4) and these 
are shown in Fig. ED Also shown here are the (non-PBH) constraints associated with the spectral dis- 
tortions in the cosmic microwave background induced by the dissipation of intermediate scale density 
perturbations and the COBE quadrupole measurement. This shows that one needs the fluctuation 
amplitude to decrease with increasing scale in order to produce PBHs and the lines corresponding to 
various slopes in the e(M) relationship are also shown in Fig. |2~1 

4-2. Refinements of Simplistic Analysis 

The criterion for PBH formation given above is rather simplistic and needs to be tested with 
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by Nadezhin et al. |136j . These roughly confirmed the criterion 5 > 7 for PBH formation, although 
the PBHs could be be somewhat smaller than the horizon. In recent years several groups have carried 
out more detailed hydrodynamical calculations and these have refined the 6 > 7 criterion, suggesting 
that one needs 5 > 0.7 rather than 5 > 0.3 |139| I157j . This in turn affects the estimate for (3{M) 
given by eqn (4). 

A particularly interesting development has been the application of "critical phenomena" to 
PBH formation. Studies of the collapse of various types of spherically symmetric matter fields have 
shown that there is always a critical solution which separates those configurations which form a black 
hole from those which disperse to an asymptotically flat state. The configurations are described by 
some index p and, as the critical index p c is approached, the black hole mass is found to scale as 
(p ~ Pc) v for some exponent rj. This effect was first discovered for scalar fields [37| but subsequently 
demonstrated for radiation jHSj and then more general fluids with equation of state p = 7/? |129[ ll()7j . 

In all these studies the spacetime was assumed to be asymptotically flat. However, Niemeyer & 
Jedamzik |138j have applied the same idea to study black hole formation in asymptotically Friedmann 
models and have found similar results. For a variety of initial density perturbation profiles, they find 
that the relationship between the PBH mass and the the horizon-scale density perturbation has the 
form 

M = KM H (6 - 8 c y (8) 

where Mh is the horizon mass and the constants are in the range 0.34 < 77 < 0.37, 2.4 < K < 11.9 
and 0.67 < 5 C < 0.71 for the various configurations. More recently, Musco et al. |135j have found that 
the critical overdensity lies in the lower range 0.43 < 5 C < 0.47 if one only allows growing modes at 
decoupling (which is more plausible if the fluctuations derive from inflation). They also find that the 
exponent 77 is modified if there is a cosmological constant. Since M — > as 5 — > S c , the existence of 
critical phenomena suggests that PBHs may be much smaller than the particle horizon at formation 
and this also modifies the mass spectrum |168| 1103"! I67| IB^j. Although Hawke & Stewart [TH] claim 
that the formation of shocks prevents black holes forming on scales below 10~ 4 of the horizon mass, 
this has been disputed |135j . 

It should be stressed that the description of fluctuations beyond the horizon is somewhat 
problematic and it is clearer to use a gauge-invariant description which involves the total energy or 
metric perturbation |157j . Also the derivation of the mass spectrum given by eqn (4) is based on 
Press-Schechter theory and it is more satisfactory to use peaks theory. Both these points have been 
considered by Green et al. [70]. They find that that the critical value for the density contrast is 
around 0.3, which is close to the value originally advocated 30 years ago! 

Another refinement of the simplistic analysis which underlies eqn (4) concerns the assumption 
that the fluctuations have a Gaussian distribution. Bullock & Primack [22] and Ivanov jHH] have 
pointed out that this may not apply if the fluctuations derive from an inflationary period (discussed 
in Section 4.3). So long as the fluctuations are small (50/0 <C 1), as certainly applies on a galactic 
scale, this assumption is valid. However, for PBH formation one requires 50/0 ~ 1, and, in this case, 
the coupling of different Fourier modes destroys the Gaussianity. Their analysis suggests that ft(M) 
can be very different from the value indicated by eqn (4) but it still depends very sensitively on e. 

4-3. PBHs and Inflation 

Inflation has two important consequences for PBHs. On the one hand, any PBHs formed 
before the end of inflation will be diluted to a negligible density. Inflation thus imposes a lower limit 
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on the PBH mass spectrum: 

M > M min = M P (T RH /T P )- 2 (9) 

where Trh is the reheat temperature and Tp m 10 19 GeV is the Planck temperature. The CMB 
quadrupole measurement implies Trh ~ 10 16 GeV, so M min certainly exceeds 1 g. On the other hand, 
inflation will itself generate fluctuations and these may suffice to produce PBHs after reheating. If 
the inflaton potential is V ((/)), then the horizon-scale fluctuations for a mass-scale M are 

^(51 (io) 

where a prime denotes d/d<j) and the right-hand side is evaluated for the value of when the mass-scale 
M falls within the horizon. In the standard chaotic inflationary scenario, one makes the "slow-roll" 
and "friction-dominated" assumptions: 

£ = {M P V'/Vf < 1, 77 = MpV"/V < 1. (11) 

Usually the exponent n characterizing the power spectrum of the fluctuations, \Sk\ 2 ~ k n , is very close 
to but slightly below 1: 

n=l + 4£- 277 «l. (12) 

omce e scales as M^ 1- ™^ 4 , this means that the fluctuations are slightly increasing with scale. The 
normalization required to explain galaxy formation (e m 10~ 5 ) would then preclude the formation of 
PBHs on a smaller scale. If PBH formation is to occur, one needs the fluctuations to decrease with 
increasing mass (n > 1) and, from eqn (12), this is only possible if the scalar field is accelerating 
sufficiently fast that [2E] 

V"/V > (1/2)(V'/V) 2 . (13) 

This condition is certainly satisfied in some scenarios jH2] and, if it is, eqn (4) implies that the 
PBH density will be dominated by the ones forming immediately after reheating. This is because the 
volume dilution of the PBHs forming shortly before the end of inflation will dominate the enhancement 
associated with eqn (4). However, it should be stressed that the validity of eqn (10) at the very end 
of inflation is questionable since the usual assumptions may fail then |1U81 1118j . 

Since each value of n corresponds to a straight line in FigE3 any particular value for the reheat 
time t\ corresponds to an upper limit on n. This limit is indicated in Fig|H]l which is taken from [3*2*] . 
apart from a correction pointed out by Green & Liddle (HHl- Similar constraints have been obtained 
by several other people jl()2| ITT] . The figure also shows how the constraint on n is strengthened if the 
reheating at the end of inflation is sufficiently slow for there to be a dust-like phase [221 EE] • However, 
it should be stressed that not all inflationary scenarios predict that the spectral index should be 
constant. Indeed, Hodges & Blumenthal [SHI have pointed out that one can get any spectrum for the 
fluctuations by suitably choosing the form of V(4>). For example, eqn (10) suggests that one can get a 
spike in the spectrum by flattening the potential over some mass range (since the fluctuation diverges 
when V goes to 0). This idea was exploited by Ivanov et al. [^U] , who fine-tuned the position of 
the spike so that it corresponds to the mass-scale associated with microlensing events observed in the 
Large Magellanic Cloud [I]. 

It should be noted that the relationship between the variance of the mass fluctuations relevant 
for PBHs and the present day horizon-scale density fluctuations in the inflation scenario is not trivial. 
For scale-free fluctuations, there is a simple relationship between them but it is more complicated 
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Fig. 3. Constraints on spectral index n in terms of reheat time t\ 



otherwise. Blais et al. fH] nn d that the mass fluctuations are reduced by 34% for a spectral index 
in the range 1 < n < 1.3, while Polarski [144] finds a further 15% reduction if there is a cosmological 
constant. 

Even if PBHs never actually formed as a result of inflation, studying them places important 
constraints on the many types of inflationary scenarios. Besides the chaotic scenario discussed above, 
there are also the variants of inflation described as designer [HEl EH 1169] , supernatural |147j , super- 
symmetric jHS], hybrid jSHUHII], multiple |167j . oscillating |162j . ghost |134j . running mass [1 14] and 
saddle [52] • There are also scenarios in which the inflaton serves as dark matter |116j . PBH formation 
has been studied in all of these models. Note that in the standard scenario inflation ends by the decay 
of the inflaton into radiation. However, in the preheating scenario inflation ends more rapidly because 
of resonant coupling between the inflaton and another scalar field. This generates extra fluctuations, 
which are not of the form indicated by eqn (10) and these might also produce PBHs |5T| I14[ 157)1 I69j. 
However, a recent analysis by Suyama et al. |161j . incorporating a full non-linear lattice simulation 
for various models, suggests that PBH formation is insignificant. There have also been a variety of 
other proposals for generating perturbations at the end of inflation |1U91 11171 1152j and the probability 
of PBH formation would then again be unrelated to eqn (10). 

5. PBHs and Dark Matter 

Roughly 30% of the total density of the Universe is now thought to be in the form of "cold 
dark matter". Recently there has been a lot of interest in whether PBHs could provide this, since 
those larger than 10 15 g would not have evaporated yet and would certainly be massive enough to be 
dynamically "cold" . It is also possible that the Planck relics of evaporated PBHs could provide the 
dark matter and this is discussed in Section 7.1. 
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5.1. Constraints on PBHs in the Galactic Halo 

There are some mass ranges in which PBHs are already excluded from providing dark matter in 
galactic halos. For example, femtolensing of gamma-ray bursts by PBHs precludes those in the mass 
range 10 17 — 10 20 g from having a critical density and microlensing of stars in the Large Magellanic 
Cloud, while allowing a tenth-critical-density at around a solar mass, excludes 10 26 — 10 34 g PBHs 
j3]. However, there are no constraints in the intermediate (sublunar) mass range 10 20 — 10 26 g [T7] . 
Note that LISA might detect the gravitational impulse induced by any nearby passing PBH [TJ 1156] . 
However, this method would not work below 10 14 g (because the effect would be hidden by the Moon) 
or above 10 20 g (because the encounters would be too rare). 

One possibility is that PBHs with a mass of around 1M could have formed at the quark- 
hadron phase transition at 10 _5 s because of a temporary softening of the equation of state then. If 
the QCD phase transition is assumed to be 1st order, then hydrodynamical calculations show that 
the value of 5 required for PBH formation is reduced below the value which pertains in the radiation 
era [HI]- This means that PBH formation will be strongly enhanced at the QCD epoch, with the 
mass distribution peaking at around the horizon mass then I166| I143j . Such PBHs might be able 
to explain the MACHO microlensing results and also the claimed microlensing of quasars [77]. 
Another possibility is that PBHs with a mass of around 1O~ 7 M could form in TeV quantum gravity 
scenarios [*"7j . 

One of the interesting implications of these scenarios is the possible existence of a halo popu- 
lation of binary black holes |137j . With a full halo of such objects, there could be a huge number of 
binaries inside 50 kpc and some of these could be coalescing due to gravitational radiation losses at 
the present epoch [!"""]. If the associated gravitational waves were detected, it would provide a unique 
probe of the halo distribution [HE]- Gravity waves from binary PBHs would be detectable down to 
1O" 5 M using VIRGO, 1O~ 7 M using EURO and 1O- U M using LISA |H7j. 

5.2. Large-Scale Structure and Supermassive Black Holes 

PBHs might also play a role in the formation of large-scale structure [4*9*] 11491 ESI HOlj • For 
example, it has been realized for some time that Poisson fluctuations in the number density of PBHs 
can generate appreciable density perturbations on large scales if the PBHs are big enough |1311 l2"t""[ 
E7J EH] • This idea has recently been invoked to explain voids and Lyman-alpha clouds [2j • In such 
scenarios, it is important to know how much a PBH can grow through accretion and it has been 
argued that this could be significant even in the period after matter-radiation equality [125 . 

Several people have suggested that the 10 6 — 1O 8 M supermassive black holes thought to 
reside in galactic nuclei could be of primordial origin. For example, a 1st order phase transition could 
produce clusters of PBHs, around which a single supermassive black hole might then condense [IB] • 
Alternatively, it has been proposed that inflation could produce closed domain walls, which then 
collapse to form 1O 8 M PBHs |l()lj ). Another scheme invokes more modest mass PBHs of 1O 3 M , 
resulting from a feature in the inflaton potential, which then grow to 1O 8 M by accretion [20]. A 
related idea comes from Bean and Magueijo, who have suggested [1*5] that PBHs may accrete from 
the quintessence field which is invoked to explain the acceleration of the Universe. However, they use 
a Newtonian formula for the accretion rate |171j and this is known to be questionable [27] . Recent 
studies of the accretion of a scalar field by a PBH also indicate that this is unlikely [73*] I74j . 
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6. Cosmic Rays from PBHS 



A black hole of mass M will emit particles like a black-body of temperature [SO] 
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GeV. 



(14) 



This assumes that the hole has no charge or angular momentum. This is a reasonable assumption 
since charge and angular momentum will also be lost through quantum emission but on a shorter 
timescale than the mass |141j . This means that it loses mass at a rate 



where the factor f{M) depends on the number of particle species which are light enough to be emitted 
by a hole of mass M, so the lifetime is 



The factor / is normalized to be 1 for holes larger than 10 17 g and such holes are only able to emit 
"massless" particles like photons, neutrinos and gravitons. Holes in the mass range 10 15 g < M < 
10 17 g are also able to emit electrons, while those in the range 10 14 g < M < 10 15 g emit muons which 
subsequently decay into electrons and neutrinos. The latter range includes, in particular, the critical 
mass M* for which r equals the age of the Universe. If one assumes the currently favoured age of 13.7 
Gyr, this corresponds to M* = 5 x 10 14 g |124j . 

Once M falls below 10 14 g, a black hole can also begin to emit hadrons. However, hadrons are 
composite particles made up of quarks held together by gluons. For temperatures exceeding the QCD 
confinement scale of Aqcd = 250 — 300 GeV, one would therefore expect these fundamental particles 
to be emitted rather than composite particles. Only pions would be light enough to be emitted below 
Aqcd- Since there are 12 quark degrees of freedom per flavour and 16 gluon degrees of freedom, 
one would also expect the emission rate (i.e. the value of /) to increase dramatically once the QCD 
temperature is reached. One can regard the black hole as emitting quark and gluon jets of the kind 
produced in collider events |12()j . The jets will decay into hadrons over a distance which is always 
much larger than the size of the hole, so gravitational effects can be neglected. The hadrons may 
then decay into astrophysically stable particles through weak and electomagnetic decays. 

To find the final spectra of stable particles emitted from a black hole, one must convolve 
the Hawking emission spectrum with the jet fragmentation function. This gives the instantaneous 
emission spectrum shown in Fig. 4 for a T = 1 GeV black hole |122j . The direct emission just 
corresponds to the small bumps on the right. All the particle spectra show a peak at 100 MeV due to 
pion decays; the electrons and neutrinos also have peaks at 1 MeV due to neutron decays. In order to 
determine the present day background spectrum of particles generated by PBH evaporations, one must 
first integrate over the lifetime of each hole of mass M and then over the PBH mass spectrum |122j . 
In doing so, one must allow for the fact that smaller holes will evaporate at an earlier cosmological 
epoch, so the particles they generate will be redshifted in energy by the present epoch. 

If the holes are uniformly distributed throughout the Universe, the background spectra should 
have the form indicated in FigED All the spectra have rather similar shapes: an E~ 3 fall-off for 
E > 100 MeV due to the final phases of evaporation at the present epoch and an E~ l tail for 
E < 100 MeV due to the fragmentation of jets produced at the present and earlier epochs. Note that 



M 



5 x 10 25 (M /g)~ 2 f{M) g s 



(15) 



r(M) = 6 x 10' 27 f(M)-\M/g) 3 s. 



(16) 
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Fig. 4. Instantaneous emission from a 1 GeV black hole 



the E" 1 tail generally masks any effect associated with the mass spectrum of smaller PBHs which 
evaporated at earlier epochs (221- The situation is more complicated if the PBHs evaporating at the 
present epoch are clustered inside our own Galactic halo (as is most likely). In this case, any charged 
particles emitted after the epoch of galaxy formation (i.e. from PBHs only somewhat smaller than 
M*) will have their flux enhanced relative to the photon spectra by a factor £ which depends upon the 
halo concentration factor and the time for which particles are trapped inside the halo by the Galactic 
magnetic field. This time is rather uncertain and also energy-dependent. At 100 MeV one has £ ~ 10 3 
for electrons or positrons and £ ~ 10 4 for protons and antiprotons. MacGibbon & Carr |121j first 
used the observed cosmic ray spectra to constrain Qpbk but their estimates have been subsequently 
refined. 




-3 -2 -i o 

Log l0 (Kinetic Energy /GeV) 



Fig. 5. Spectrum of particles from uniformly distributed PBHs 
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6.1. Gamma-Rays 

EGRET observations |158j of the 7-ray background between 30 MeV and 120 GeV leads to 
an upper limit I^pbh < (5.1 ± 1.3) x 10~ 9 h~ 2 , where h is the Hubble parameter in units of 100 [3*U] . 
This is a refinement of the original Page-Hawking limit, but the form of the spectrum suggests that 
PBHs do not provide the dominant contribution. If PBHs are clustered inside our own Galactic 
halo, then there should also be a Galactic 7-ray background and, since this would be anisotropic, 
it should be separable from the extragalactic background. The ratio of the anisotropic to isotropic 
intensity depends on the Galactic longtitude and latitude, the Galactic core radius and the halo 
flattening. Wright claims that such a halo background has been detected [165J. His detailed fit to 
the EGRET data, subtracting various other known components, requires the PBH clustering factor 
to be (2 — 12) x 10 5 /i _1 , comparable to that expected. 

6.2. Antiprotons and Antideuterons 

Since the ratio of antiprotons to protons in cosmic rays is less than 10 -4 over the energy range 
100 MeV— 10 GeV, whereas PBHs should produce them in equal numbers, PBHs could only contribute 
appreciably to the antiprotons j!64j . It is usually assumed that the observed antiproton cosmic rays 
are secondary particles, produced by spallation of the interstellar medium by primary cosmic rays. 
However, the spectrum of secondary antiprotons should show a steep cut-off at kinetic energies below 
2 GeV, whereas the spectrum of PBH antiprotons should increase with decreasing energy down to 
0.2 GeV. Also the antiproton fraction should tend to 0.5 at low energies, so these features provide 
a distinctive signature |103j . MacGibbon & Carr originally calculated the PBH density required to 
explain the interstellar antiproton flux at 1 GeV and found a value somewhat larger than the 7-ray 
limit |121j . More recent data on the antiproton flux below 0.5 GeV come from the BESS balloon 
experiment |170j and Maki et al. |128j tried to fit this in the PBH scenario by using Monte Carlo 
simulations of cosmic ray propagation. A more recent attempt to fit the antiproton data comes from 
Barrau et al. [H] and is shown in FigHT]PBHs might also be detected by their antideuteron flux 



PBH antiproton spectra 




antiproton kinetic energy (GeV) 



Fig. 6. Comparison of PBH emission and antiproton data from Barrau et al. 
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6. 3. PBH Explosions 

One of the most striking observational consequences of PBH evaporations would be their final 
explosive phase. However, in the standard particle physics picture, where the number of elementary 
particle species never exceeds around 100, the likelihood of detecting such explosions is very low. 
Indeed, in this case, observations only place an upper limit on the explosion rate of 5 x 10 8 pc _3 y _1 
[""l ll51j . This compares to Wright's 7-ray halo limit of 0.3 pc _3 y _1 [165j and the Maki et al. antiproton 
limit of 0.02 pc^y" 1 

However, the physics at the QCD phase transition is still uncertain and the prospects of 
detecting explosions would be improved in less conventional particle physics models. For example, in 
a Hagedorn-type picture, where the number of particle species exponentiates at the the quark-hadron 
temperature, the upper limit is reduced to 0.05 pc~ 3 y _1 jSHJ - Cline and colleagues have argued that 
one might expect the formation of a QCD fireball at this temperature EH] and this might even 
explain some of the short period 7-ray bursts observed by BATSE. They claim to have found 42 
candidates of this kind and the fact that their distribution matches the spiral arms suggests that 
they are Galactic [H]. Although this proposal is speculative and has been disputed jHHj, it has the 
attraction of making testable predictions (eg. the hardness ratio of the burst should increase as 
the duration decreases). A rather different way of producing a 7-ray burst is to assume that the 
outgoing charged particles form a plasma due to turbulent magnetic field effects at sufficiently high 
temperatures 

Some people have emphasized the possibility of detecting very high energy cosmic rays from 
PBHs using air shower techniques [72*1 l4"4t flllj . However, this is refuted by the claim of Heckler [HI] 
that QED interactions could produce an optically thick photosphere once the black hole temperature 
exceeds T crit = 45 GeV. In this case, the mean photon energy is reduced to m e (Te# /T cr it) 1 ^ , which 
is well below Tbh, so the number of high energy photons is much reduced. He has proposed that a 
similar effect may operate at even lower temperatures due to QCD effects [""I"] . Several groups have 
examined the implications of this proposal for PBH emission [521 ETJ- However, these arguments 
should not be regarded as definitive since MacGibbon et al. claim that QED and QCD interactions 
are never important 



7. PBHs as a Probe of Quantum Gravity 

7.1. Planck Mass Relics 

Some people have speculated that black hole evaporation could cease once the hole gets close to 
the Planck mass |12(i[ "2T" [ "4~*"]. For example, in the standard Kaluza-Klein picture, extra dimensions are 
assumed to be compactified on the scale of the Planck length. This means that the influence of these 
extra dimensions becomes important at the energy scale of 10 19 GeV. Such effects could conceivably 
result in evaporation ceasing at the Planck mass. Various non-quantum-gravitational effects (such 
as higher order corrections to the gravitational Lagrangian or string effects) could also lead to stable 
relics j""2] but the relic mass is usually close to the Planck scale. 

Another possibility, as argued by Chen & Adler [""""j] > i s that stable relics could arise if one 
invokes a "generalized uncertainty principle". This replaces the usual uncertainty principle with one 
of the form 

**>l p+ llf, (IT) 
where the second term is supposed to account for self-gravity effects. This means that the black hole 
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temperature becomes 

This reduces to the standard Hawking form for M ^> Mp but it remains finite instead of diverging at 
the Planck mass itself. Cavaglia and colleagues have extended this argument to higher dimensional 
black holes [22] and also included the effects of thermal fluctuations [3*lj . 

Whatever the cause of their stability, Planck mass relics would provide a possible cold dark 
matter candidate |119j . Indeed this leads to the "relics" constraints indicated in Fig.l, Fig. 2 and 
Fig.3. In particular, such relics could derive from inflationary PBHs [T3*l ll7|. If the relics have a mass 
nMp, then the requirement that they have less than the critical density implies [22] 

(3{M) < l$~ 27 K-\M/M P f/ 2 (19) 

for the mass range 

(T RH /T P y 2 M Pl < M < 10 u k 2/5 M p . (20) 

The upper mass limit arises because PBHs larger than this dominate the total density before they 
evaporate. Producing a critical density of relics obviously requires fine-tuning of the index n. Also 
one needs n wl.3, which is barely compatible with the WMAP results Nevertheless, this is 

possible in principle. In particular, it has been argued that hybrid inflation could produce relics from 
10 5 g PBHs formed at l(T 32 s [HHllTfiS]. 



7.2. PBHs and Brane Cosmology 

In "brane" versions of Kaluza-Klein theory, some of the extra dimensions can be much larger 
than the usual Planck length and this means that quantum gravity effects may become important at 
a much smaller energy scale than usual. If the internal space has n dimensions and a compact volume 
V n , then Newton's constant Gn is related to the higher dimensional gravitational constant Go and 
the value of the modified Planck mass Mp is related to the usual 4-dimensional Planck mass M 4 by 
the order-of-magnitude equations: 

G N ~ G D /V n , M n P +2 ~ M 2 jV n . (21) 

The same relationship applies if one has an infinite extra dimension but with a "warped" geometry, 
provided one interprets V n as the "warped volume". In the standard model, V n ~ I/M4 and so 
Mp ~ M4. However, with large extra dimensions, one has V n ^> 1/M% and so Mp <C M4. 

Brane cosmology modifies the standard PBH formation and evaporation scenario. In partic- 
ular, in Randall- Sundrum type II scenarios, a p 2 term appears in the Friedmann equation and this 
dominates the usual p term for 

T > 10 18 (Z/Z P r 1/4 GeV (22) 

where I is the scale of the extra dimension. This exceeds 1 TeV for I < 0.2 mm. Black holes which 
are smaller than I are effectively 5-dimensional. This means that they are cooler and evaporate more 
slowly than in the standard scenario: 

T BH oc M" 1/2 , r oc M 2 . (23) 
The critical mass for PBHs evaporating at the present epoch is now in the range 3 x 10 9 g to 4 x 10 14 g. 

This mnHifipR thp standard PBH pvannratinn rnnstraintfi 171 1 13811 and if snmp rnsmir ravs pnmp frnm 
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PBHs, it means that these could probe the extra dimensions [1541 1155] . Another complication is that 
accretion could dominate evaporation during the era when the p 2 term dominates the p term in the 
Friedmann equation |13Uj . 

7. 3. Black Hole Production at Accelerators 

One exciting implication of these scenarios is that quantum gravitational effects may arise at 
the experimentally observable TeV scale. If so, this would have profound implications for black hole 
formation and evaporation since black holes could be generated in accelerator experiments, such as 
the Large Hadron Collider (LHC) |4*7l IBTl Il50| . Although the black holes produced in accelerators 
should not themselves be described as "primordial" , since they do not form in the early Universe, it 
is clear that these experiments will also have profound implications for PBHs because, at sufficiently 
early times, the effects of the extra dimensions will be cosmologically important (see Section 7.2). 

Two partons with centre-of-mass energy y/s will form a black hole if they come within a 
distance corresponding to the Schwarzschild radius rs for a black hole whose mass Mbh is equivalent 
to that energy jHJ- Thus the cross-section for black hole production is 

a BH « 7rr|e(Vi - Mg£) (24) 

where M^jp is the mass below which the semi-classical approximation fails. Here the Schwarzschild 
radius itself depends upon the number of internal dimensions: 

1 /M RH \ 1/{1+n) 
r s*TT- ( -77^ ) > ( 25 ) 



M P \ M P 

so that apH s 2 ^ n+l \ This means that the cross-section for black hole production in scattering 
experiments goes well above the cross-section for the standard model above a certain energy scale 
and in a way which depends on the number of extra dimensions. 

The evaporation of the black holes produced in this way will produce a characteristic signature 
|4"7|I61|IT5U] because the temperature and lifetime of the black holes depend on the number of internal 
dimensions: 

n + l 1 /M B/ A (n+3)/(n+1) 



TflH ~-^-' TBH ~Mp-\-Mp-) ■ (26) 

Thus the temperature is decreased relative to the standard 4-dimensional case and the lifetime is 
increased. The important qualitative effect is that a large fraction of the beam energy is converted 
into transverse energy, leading to large-multiplicity events with many more hard jets and leptons than 
would otherwise be expected. In principle, the formation and evaporation of black holes might be 
observed by the LHC when it turns on in 2007 and this might also allow one to experimentally probe 
the number of extra dimensions. On the other hand, this would also mean that scattering processes 
above the Planck scale could not be probed directly because they would be hidden behind a black 
hole event horizon [60J. 

Similar effects could be evident in the interaction between high energy cosmic rays and at- 
mospheric nucleons. Nearly horizontal cosmic ray neutrinos would lead to the production of black 
holes, whose decays could generate deeply penetrating showers with an electromagnetic component 
substantially larger than that expected with conventional neutrino interactions. Several authors have 
studied this in the context of the Pierre Auger experiment, with event rates in excess of one per 
year being predicted E3J 1150] . Indeed there is a small window of opportunity in which Auger 

miffht rlfitp.nt surh pvpnts hpfnrp thp LHC! On thp nthpr hand it must pantinnpH that thpsp pstimatps 
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may be overly optimistic. Others have claimed that black hole production in accelerators may be 
unobservable even if TeV quantum gravity is correct j3] . 

8. Conclusions 

We have seen that PBHs could be used to study the early Universe, gravitational collapse, 
high energy physics and quantum gravity. In the "early Universe" context, useful constraints can 
be placed on inflationary scenarios. In the "gravitational collapse" context, the existence of PBHs 
could provide a test of critical phenomena. In the "high energy physics" context, information may 
come from cosmic ray and even gamma-ray bursts if suitable physics is invoked at the QCD phase 
transition. In the "quantum gravity" context, the formation and evaporation of small black holes 
could be observable in cosmic ray events and accelerator experiments if the quantum gravity scale is 
around a TeV. 
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